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We systematically derive an effective Hamiltonian for the dipolar magnet LiYi_ x Ho x F4, including 
quantum corrections which arise from the transverse dipolar and hyperfine interactions. These 
corrections are derived using a generalised Schrieffer- Wolff transformation to leading order in the 
small parameters given by the ratio of the interaction energies to the energy of the first excited 
electronic state of the Holmium ions. The resulting low-energy Hamiltonian involves two-level 
systems, corresponding to the low-lying electronic states of the Holmiums, which are coupled to one 
another and to the Holmium nuclei. It differs from that obtained by treating the electronic states 
of the Holmium as a spin-1/2 with an anisotropic g-factor. It includes effective on-site transverse 
fields, and both pairwise and three-body interactions among the dipoles and nuclei. We explain the 
origins of the terms, and give numerical values for their strengths. 
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I. INTRODUCTION 

The rare-earth compound LiYi_ x Ho x F4has often been 
described as a model quantum magnet, and as such 
has been studied for over three decades^ 2 - 3 -^^. The 
basic model for this material is a diluted Ising model 
with dipole interactions, and depending on the dilu- 
tion the low-temperature phase is expected to be ei- 
ther a ferromagnet or a spin-glassi Applying a strong 
transverse field to the material introduces quantum 
fluctuations which can lead to domain wall tunnelling 
in the ferromagnet^, and to quantum melting of the 
ferromagnet^ and the spin-glass^ i 10 ' 11 ^ 2 ' 13 . Even in 
the absence of an applied transverse field, however, 
quantum tunnelling of magnetic dipoles can be ob- 
served in LiYi_ x Ho x F^ii^. Furthermore, samples with 
a Holmium concentration of 4.5% do not appear to be- 
have as a spin-glass on cooling, as expected for a clas- 
sical Ising model. Instead they show unexplained "anti- 
glass" properties which have been attributed to quan- 
tum mechanics^ This attribution is supported by the- 
oretical predictions of the static susceptibility, which 
agree with experiments only once quantum corrections 
are included^. 

Although quantum effects are implicated in several 
phenomena in LiYi_ x Ho x F4in the absence of an applied 
transverse field, there is no derivation in the literature 
of an effective low-energy Hamiltonian which contains 
them. Here we derive such a Hamiltonian, projecting 
out the high-energy electronic states to obtain a theory 
describing the low-lying electronic doublet and the nuclei. 
We shall find a Hamiltonian which has a different form 
from that appropriate in a strong applied fielcU^, and 
whose electronic part differs from the two-level model 
previously proposed for the zero-field case^. 



II. MODEL 

The magnetic degrees of freedom in LiYi_ x Ho x F4are 
the / electrons on the Ho 3+ ions. The strong spin-orbit 
coupling of the ions leads to a well-defined J = 8 for 



the ions, and an associated dipole moment \i = gLfJ-B, 
with the Lande g-factor g L = 5/4. The 2 J + 1 = 17-fold 
degeneracy of the free ion is broken by the crystal-field 
Hamiltonian, leaving a degenerate ground-state doublet. 
All matrix elements of J x and J v are zero within this 
ground-state subspace, but there are non-zero matrix el- 
ements for J z . This is the source of the strong Ising 
anisotropy in the interactions and response to an applied 
field. 

The interactions between Holmium ions are dipolar, 
giving an interaction Hamiltonian 
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The crudest way to obtain a low-energy effective Hamil- 
tonian from ((T|) is to truncate the electronic state-space 
to the ground-state doublet on each Holmium ion. It is 
possible to choose a basis for the doublet in which J 2 has 
no off-diagonal matrix elements, while the diagonal ma- 
trix elements are a and —a. With this choice of basis, 
truncating leads to the dipolar Ising model 



E 

hi 



(l-3(f v ,k) 2 Ka? (2) 



where a x ,a v and a z are the usual Pauli matrices. As 
discussed in previous studies of this system 3 - ^ 10 ' 15 the 
contact hyperfine interaction between the Holmium nu- 
clei and the electronic degrees of freedom also plays an 
important role, as it is typically of similar strength to 
the dipolar interaction. Under the simple two-state trun- 
cation scheme described above the hyperfine interaction 
also takes a simple Ising form, 



H hyp = J^AjIi-Ji (3) 
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as do other couplings to the Holmium moments. 
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The procedure of projecting into the low-energy sub- 
space of the single-ion Hamiltonian has been applied in 
the presence of an applied transverse field.— The field 
changes the low-energy states, leading to finite matrix 
elements for J x and J y within the low-energy doublet. 



This leads to corrections to the Ising forms (|2|4p that are 
pair interactions involving at least one of a x ,a v ,I x ,I v , 
and to effective field terms. 

To obtain quantum terms in the absence of an applied 
transverse field we must go beyond a simple projection 
onto the low-lying single-ion states. We thus now con- 
sider the lowest three levels of the crystal-field Hamil- 
tonian, denoting the two states of the doublet as | f), 

I) and the first excited state, with energy A, as \T). 
The parameters in such a three-level model are A and 
the matrix elements of the angular momentum operators, 
and can be obtained from a numerical diagonalisation of 
the full crystal field Hamiltonian. With an appropriate 
choice of basis this gives,— 
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(T \ Jx\T) = 2.4 = p, (6) 

-(UJy\r) = i P , (7) 
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All other matrix elements of angular momentum among 
the three states vanish. 



III. DERIVATION OF THE EFFECTIVE 
HAMILTONIAN 



Quantum corrections to the Ising Hamiltonian (|2|4j) 
appear once one includes the state at |T) because 
the interaction terms ([I]) and ([3]) couple |T) to the 
electronic ground-state doublet. The typical nearest- 
neighbour dipolar interaction energy scale is ~ 300mK 
and Aj m 38mK^ii£. Since these interaction scales are 
small compared with A the couplings to the |T) states 
can be treated in perturbation theory, and an effective 
Ising model obtained using the standard Schrieffer- Wolff 
procedure^. The resulting model will include quan- 
tum terms arising from virtual transitions between the 
ground-state doublet and the |T) state. 

To derive the effective Hamiltonian we begin by writing 
the three- level model as H = Hq + Ht, where Hq con- 
tains terms which do not couple the doublet to and 
Ht contains those which do. Measuring energy from the 
electronic ground-state doublet, Hq contains the crystal 
field term V c = A^V ) {r^ | and the Ising parts of the 
interactions, while Ht contains the parts of the interac- 
tions which involve J x and J y . Including the hyperfme 
interaction the Hamiltonian can be spilt as 
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where v, /x = x, y, z, but the prime on the sum indicates 
that we must exclude all terms that have v = fi = z. 

We then seek a unitary transformation H — * e s He~ s 
which decouples the electronic doublet from the |T) state 
to first order in Ht- Such a transformation obeys 



[S, Hq] — — Ht- 
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An S which approximately satisfies (fTTj) can be con- 
structed from Ht using projection operators, 
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The form (p~2|) for S actually obeys the relation [S, V c ) = 
—Ht, and although we have successfully eliminated the 
linear term in Ht the remainder of the commutator 
\S, H — V c ] generates new couplings between the doublet 
and |r). However these tunnelling terms are of order 
(interaction) 2 /A, and as they can only couple between 
the ground states and |r), they can only contribute to 
the effective low-energy Hamiltonian in higher order per- 
turbation theory. Therefore they do not contribute to 
the effective low-energy Hamiltonian to leading order in 
(interaction) 2 /A, and vanish when we project onto the 
Ising basis at the end of the Schrieffer- Wolff procedure. 

Having eliminated Ht, the effective two-state Hamil- 
tonian is given to lowest order in (interactions) 2 /A by^ 
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This form extends the Ising interaction Hamiltonian by 
including second-order processes in which Ht causes vir- 
tual transitions from the electronic doublet to |T) and 
then back again. 



IV. EFFECTIVE HAMILTONIAN 

Substituting the form (fT2"|) for S into Eq. [TH] and dis- 
carding irrelevant energy shifts we find that the effective 
Hamiltonian takes the form 
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This Hamiltonian operates in the space with sixteen 
states for each ion: two low-lying electronic states, each 
with the eight nuclear states. The various contributions 
to H c ff are labelled according to the interactions that give 
rise to them. 



A. Dipolar processes 

The first term in (fT7|) arises from the second-order 
dipolar processes that involve only single pairs of spins, 
and is 

Hu = »>; + h{aV) + E A^<af (18) 

We see that the pairwise dipolar interaction thus gener- 
ates both transverse interaction terms and effective mag- 
netic fields. In appendix [5] we give expressions for the 
strengths of the parameters hi and A"'* in terms of the 
parameters of the underlying three-level model. Using 
these relations we estimate the characteristic magnitude 
of these corrections by calculating the contribution from 
a single nearest-neighbour. This yields h x « 3mK, h y w 
and A 1 ^ « 0.06mK, which may be compared to the dipo- 
lar Ising energy Jff N ~ 300mK, and the fsing hyperfine 
splitting aAj ps 210mK. Note that terms such as (fT8|) de- 
scribe tunnelling between the low-lying electronic states. 
It is important to stress that this is not sufficient to gen- 
erate tunnelling between the two electro-nuclear ground 
states of i?cr^- 

In addition to the dipolar processes involving only a 
single pair of Holmium ions there are processes involving 
three ions. In such a process the interaction between one 
pair of spins, say virtually excites the jth spin into 
r), and the dipolar interaction of this spin with the fcth 
spin brings the jth spin back into the Ising doublet. This 
generates two- and three- body interaction terms in the 
effective theory, 

Hutb = -^f E (^*Jg* + «*W ( 19 ) 

For an isolated group of three spins in which {i,j} and 
{i, k} are nearest neighbours we obtain 3 mK for the 
magnitude of the interaction terms cr^cr| , a^ai-af and 

B. Hyperfine processes 

The contributions to the effective Hamiltonian from 
the hyperfine interactions are simpler to deal with, as 



they arc confined to each site. They are 
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We see that in second order perturbation theory the 
hyperfine interaction leads to spin flipping terms which 
require the z component of the nuclear spin to change 
by two units, gives a slight correction to the longitudi- 
nal hyperfine interaction, and introduces a weak nuclear 
anisotropy. The factors outside the sum give 0.86mK for 
the interaction energy. 

C. Mixed hyperfine-dipolar processes 

Finally, the term in the effective Hamiltonian labelled 
i^ND arises due to second-order processes whereby the 
electrons are virtually excited into the state \T) by the 
dipolar interaction and then de-excited by the hyperfine 
interaction or vice versa: 

H ND = -^l^2j?f(If(l + af)+Ifay)a] 

Note that the nuclei mediate a coupling between elec- 
tronic states of the form a x ' v a z which does not occur in 

V. DISCUSSION 

In recent paper o 10 ' 11 ^ 2 ' 13 it was shown that the non- 
Ising parts of the dipolar coupling provide a route by 
which an applied transverse field can destroy spin-glass 
order in LiYi_ x Ho x F4. Although our model contains 
spontaneous transverse fields we do not expect them to 
destroy the spin-glass phase in this way, because the 
model retains time-reversal symmetry. An interesting 
feature of the effective Hamiltonian is that the spon- 
taneous transverse fields and the two- and three- body 
interactions are all correlated. Theoretical work on spin 
models with several forms of random couplings, such as 
the random transverse-field Ising model^, generally con- 
siders different interaction terms as independent. 

The approach given here could be straightfor- 
wardly extended to derive an effective Hamiltonian for 
LiYi_ x Ho x F4in an applied transverse field, so long as the 
field strength is sufficiently small for perturbation theory 
to apply. This generates a variety of new processes at 
second order. The straightforward process, in which the 
applied field excites to the |T) state and back again, leads 
to an electronic tunnelling term oc B 2 /A. But note that 
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there will also be mixed processes, for example where the 
applied field causes the transition to |T) and the dipole 
interaction causes the transition back into the electronic 
doublet, which will generate electronic tunnelling terms 
oc B/A. 

We stress that although our Hamiltonian contains 
terms which generate transitions between the states of 
the electronic Ising doublet, and further such terms will 
be introduced by an applied field, these terms alone do 
not couple the doubly-degenerate electro-nuclear ground 
states of the Ising single-ion Hamiltonian. To generate 
quantum fluctuations between the electro-nuclear ground 
states requires processes which flip the nuclear spin, and 
hence the splittings of the electro-nuclear ground states 
will be much smaller than the tunnelling terms connect- 
ing the bare electronic states. 

Although we do not expect the quantum corrections 
to formally destroy the ordered phase unless there is 
an applied field, they can still affect the thermodynam- 
ics, changing the susceptibilities and moving the phase 
boundary. However, we only expect this to occur at 
very low temperatures. In this context we note that the 
experimental susceptibility of LiYi_ x Ho x F4in the 10 - 
100 mK temperature regime of the antiglass experiment 
has been reproduced by a quantum theory of a two-level 
model^. That model is obtained by neglecting the nu- 
clei and treating the Holmium ions as spin-1/2 ions with 
an anisotropic g-factor, i.e. writing Jj~ = .g'Vf, with 
H = x,y, z. We note that the resulting Hamiltonian 
is very different from that which would be obtained by 
neglecting the nuclei in our model: there are no spon- 
taneous field terms and the non-Ising interactions de- 
cay as 1/r 3 (here l/r 6 ), while the quoted values for g 
give the energy scale for the largest non-Ising coupling of 
f=a 30mK, whereas here we have 3mK. 

While we do not expect the small corrections derived 
here to affect thermodynamics except at very low temper- 
atures, they may be relevant to understanding dynam- 
ics at much higher temperatures. Quantum tunnelling 
of the magnetisation has been observed in the magneti- 
sation relaxation and susceptibility of the dilute com- 
pound LiYo.998Hoo.oo2F4, due to both single-ion^ and 
two-ion processes^, and it would be interesting to com- 
pare the details of these results with the processes given 
here. An understanding of the single or few-ion electro- 
nuclear dynamics, based on the Hamiltonian (|17p . may 
also help to explain the antiglass experiments^. In the 
high-temperature regime these experiments show a char- 
acteristic relaxation time for the magnetisation which is 
activated, with a barrier similar to the width of the hy- 
pcrfinc multiplct and the nearest-neighbour interaction 
strength. Extrapolating this behaviour indicates that 
this particular activated dynamics freezes out on the ex- 
perimental frequency scale as the temperature is lowered 
into the antiglass regime. The quantum terms derived 
here, although small, may perhaps then provide a route 
to the observed dynamics. The potential significance of 
quantum electro-nuclear dynamics for the antiglass ex- 



periment has been noted by Atsarkini^, who proposes a 
specific relaxation mechanism due to interactions of the 
form (HI]). 

VI. CONCLUSIONS 

In this paper we have motivated the need to go beyond 
the simple dipolar electro-nuclear Ising Hamiltonian com- 
monly used to describe LiYi_ x Ho x F4, and have derived 
an effective low energy Hamiltonian which includes the 
quantum corrections caused by the transverse elements 
of the dipolar and hypcrfinc interactions. We have given 
estimates for the typical magnitudes of these correction 
terms and have shown that they are typically about one 
percent as strong as the energy scale associated with the 
Ising interactions. As a result, we do not expect any 
qualitative changes to the Ising phase diagram, but as 
we have highlighted in this paper, these quantum cor- 
rections can describe a large variety of single and many- 
body processes which might play significant roles in the 
observed dynamics. Thus, for low temperatures, our ef- 
fective Hamiltonian should serve as a good starting point 
for a microscopic investigation of the dynamical physics 
of LiYi_ x Ho x F4, and as it can be used across the whole 
dilution series, should contain the rich low energy physics 
which characterises the spin glass, free ion, and presum- 
ably, the "antiglass" phases of this material. 
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APPENDIX A: RELATION OF MICROSCOPIC 

PARAMETERS TO THE EFFECTIVE 
INTERACTION AND ON-SITE TERMS IN H D 

The components of the effective magnetic field h,, de- 
fined in Eq. [THJ are related to the original microscopic 
Hamiltonian by, 



h < = V^ ((J ^ )2 ~ (J ^ )2) 
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The components of the effective magnetic field vanish 
in the undiluted crystal, as expected from the crystal 
symmetry. 

The transverse dipolar interactions between spins are 
described by the following couplings, 
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These are the only couplings generated; there are no 
terms such as a x a z , since these would break time- reversal 
symmetry. 
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